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Introduction 


1.1.  Background 

The  multiplex  and  throughput  advantages  of  Fourier  spectroscopy  have  made 
it  possible  to  obtain  rapid  measurements  of  local  emission  spectra.  The  high 
resolution  emission  spectra  that  are  available  and  will  be  coming  from  the 
Stratospheric  Cryogenic  Interferometer  Balloon  Experiment  (SCRIBE)  (Murcray  et 
al . ,  198A,  1985)  promise  to  provide  important  atmospheric  information  un¬ 
attainable  from  solar  occultation  absorption  spectra  or  other  means.  Of 
special  interest  are  the  diurnal  variation  of  ozone  -  related  trace  gases  in  the 
stratosphere  and  anomalous  variation  with  height  of  spectra  in  band  passes 
used  for  temperature  and  compositional  sounding,  such  as  CO2  Q-branches 
subject  to  line  mixing.  The  SCRIBE  data  will  also  provide  useful  information 
for  establishing  the  background  radiation  levels  with  which  Air  Force  electro- 
optical  systems  must  contend. 

We  proposed  to  begin,  during  the  first  phase  of  interpretation  of  SCRIBE 
spectra,  three  complementary  investigations,  namely,  the  validation  of  the 
spectra,  the  retrieval  of  trace  gas  profiles,  and  the  retrieval  of  temperature 
profiles.  Contracts  were  awarded  both  to  us  and  to  OPtiMetrics,  Inc.  This  is 
fortunate,  since  there  is  much  to  be  done  with  the  spectra,  and  our  program  is 
ambitious.  Efficiency  required,  however,  that  redundancy  and  duplication  of 
effort  be  minimized.  Our  AFGL  monitors  suggested  that  OPtiMetrics  have  the 
primary  responsibility  for  calibration  and  validation,  while  we  concentrate  on 
retrievals,  particularly  of  temperature.  We  were  happy  to  accept  this  plan  of 
action . 

Prior  to  commencement  of  the  contract  research,  Professor  Aaron  Goldman 
of  the  University  of  Denver  sent  us  samples  of  observed  and  synthetic  data. 
From  an  examination  of  observed  and  synthetic  spectra  near  the  15|jm  CC>2 
Q-branch,  we  noticed  that  there  appeared  to  be  evidence  for  a  substantial  line 
coupling  effect.  The  sign  of  the  discrepancy  seemed  to  be  compatible  with 
the  departure  of  laboratory  measurements  at  room  temperature  from  synthetic 
laboratory  spectra  calculated  by  Dr.  Michael  Hoke,  of  AFGL's  Optical  Phvsics 
Division,  for  the  same  temperature,  with  and  without  line  coupling. 


-1- 


Since  this  region  of  the  spectrum  is  so  important  for  remote  sensing,  and 
since  there  were  no  observed  or  calculated  line  coupling  data  at  temperatures 
lower  than  room  temperature,  it  was  agreed  that  line  coupling  studies  using 
SCRIBE  data  would  have  high  priority  in  this  first  phase  of  the  research, 
together  with  retrieval  development  and,  of  course,  the  adaptation  of  FASC0D2 
(Clough  et  al,  1986)  for  calculation  of  synthetic  SCRIBE  spectra,  which  is 
necessary  to  carry  out  the  other  two  objectives. 


1 ■ 2  Technical  Objectives  of  Study 


This  research  proposed  investigation  of  the  potential  of  the  unique 
SCRIBE  high  resolution  spectra  to  provide  information  on  atmospheric 
temperature  and  trace  gas  concentrations  and  to  verify  and  validate 
spectroscopic  modeling  approaches  such  as  that  for  line  coupling.  It  was 
hoped  that  eventually  these  retrieval  and  spectroscopic  studies  could  be 
extended  to  observations  of  the  diurnal  variations  of  trace  gases  and 
comparison  of  SCRIBE  retrieved  concentrations  to  those  obtained  from  state-of- 
the-art  photochemical  models. 


The  technical  objectives  for  this  first  phase  of  the  study  were  thus 
threefold,  as  follows:  (1)  development  of  a  FASC0D2  routine  for  the  forward 
problem  of  calculating  simulated  radiances  received  at  the  SCRIBE  altitude  and 
angle-of -view;  (2)  using  this  routine  to  develop  a  FASC0D2 -based  temperature 
retrieval  for  SCRIBE  data,  and  (3)  investigating  the  effect  of  line  coupling 
near  the  15  pm  Q-branch  on  FASC0D2  simulated  radiance  spectra  and  verifying 
this  behavior  in  the  SCRIBE  data. 


1 . 3  Report  Overview 


Chapter  2  below  contains  a  brief  description  of  the  SCRIBE  data,  how  it 
was  obtained,  calibration  and  noise  problems,  and  the  data  subset  that  was 
used  in  our  study. 


Chapter  3  describes  the  SCRIBE  temperature  retrieval  studies,  starting 
with  the  forward  problem  calculations  and  development  of  a  FASC0D2 -based 
retrieval  scheme  adapted  for  SCRIBE.  Channel  selection  and  sensitivity 
studies  with  these  provisional  channels  are  described,  and  the  chapter  is 
concluded  with  an  examination  and  discussion  of  preliminary  retrieval  results. 


Chapter  4  discusses  the  investigation  of  the  effects  of  line  coupling  on 
SCRIBE  spectra  in  the  vicinity  of  15/um,  as  uncovered  by  our  investigation. 

The  effects  are  striking  and  may  have  important  implications  for  remote 
sensing . 

Finally,  Chapter  5  contains  a  summary,  and  recommendations  for  future 
research . 

2.  SCRIBE  Data  Set 

The  SCRIBE  instrument  is  a  cryogenic  Michelson  interferometer  spectro¬ 
meter  carried  aloft  by  balloon  for  the  purpose  of  observing  the  atmospheric 
emission  spectra.  It  was  flown  in  1983,  1984  and  1986,  each  time  obtaining 
spectra  from  600  to  1500  cm"^  with  a  spectral  resolution,  after  apodization, 
of  0.06  to  0.07  cm’^.  Each  succeeding  flight  resulted  in  spectra  with  higher 
signal-to-noise  ratio  and  easier  to  calibrate,  but  the  1986  data  were  not 
reduced  and  calibrated  in  time  for  us  to  work  with  them.  The  data  for  the 
flight  on  July  5,  1984,  which  was  launched  from  Roswell,  New  Mexico,  was 
available,  however,  and  OPtiMetrics,  Inc.,  provided  us  with  magnetic  tapes 
containing  four  "corrected"  spectra  from  float  elevation  30.5  km,  at  zenith 
angles  180  degrees  (nadir),  88.1,  93.2  and  93.7  degrees.  These  are  shown  in 
Figs.  1  to  4,  respectively,  with  the  spectral  region  from  600  to  1500  cm‘^ 
compressed  to  fit  each  spectrum  onto  a  single  page.  Much  of  the  spectra  are 
seen  to  be  quite  noisy,  even  after  correction,  especially  at  the  high 
frequency  end,  and  a  very  prominent  spike  appears  at  755  cm'^  on  the  88.1 
degree  spectra.  These  features  are  not  uncommon;  in  fact,  these  are 
presumably  among  the  best  spectra. 

An  even  more  serious  problem  for  interpreting  the  spectra  is  that  of 
calibration.  A  blackbody  source  was  flown,  but  it  seems  to  have  encountered 
drifts.  An  idea  of  the  magnitude  of  the  calibration  uncertainty  can  be 
obtained  by  comparing  the  four  spectra  at  the  common  feature  at  668  cm"^, 
corresponding  to  the  peak  of  the  15  ^m  CC>2  Q-branch.  Here  the  atmosphere  is 
so  opaque  at  all  angles  that  the  measured  radiances  should  be  the  same, 
corresponding  to  a  brightness  temperature  equal  to  the  ambient  temperature  at 
float  altitude,  which  was  measured  by  the  White  Sands,  New  Mexico,  radiosonde 
to  be  very  close  to  230K.  It  can  be  seen,  say  by  laying  transparenc ies  of  th 
spectra  one  on  the  other,  that  the  quas i - horizontal lv - viewed  radiances  agree 


Figure  1  Corrected  SCRIBE  data,  nadir  view,  600-1500 


Figure  3  Corrected  SCRIBE  data,  93.2  deg 
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with  each  other  (the  corresponding  brightness  temperatures  are  about  223K)  but 
the  radiance  for  nadir  view  is  much  higher,  corresponding  to  a  brightness 
temperature  of  about  250K.  This  can  be  seen  in  more  detail  in  Figs.  23  and 
30,  below.  Clearly  great  care  should  be  taken  with  the  calibration  in  inter¬ 
preting  the  spectra,  especially  if  the  angular  variation  is  used  together  with 
the  spectral  variation  in  retrievals. 

In  this  study  a  small  subset  of  the  1984  SCRIBE  data  was  used.  The 
spectral  regions  used  were  in  the  15  jjm  CO2  vibration- rotation  band,  more 
particularly  between  600  to  670  cm"^,  since  the  retrieval  studies  and  line 
coupling  investigations  fell  within  that  spectral  interval. 

3.  SCRIBE  Temperature  Retrieval  Studies 

3 . 1  Development  of  FASC0D2 - Based  Physical  Retrieval  for  SCRIBE 

Since  no  statistics  were  available  relating  radiances  at  SCRIBE  level  to 
atmospheric  state  variables,  and  since  it  would  have  required  a  prohibitive 
effort  to  acquire  a  set  of  synthetic  statistics,  a  physical  relaxation  method 
was  used  in  our  temperature  retrieval  routine.  This  required  forward 
computation  of  radiances.  The  availability  of  the  AFGL  line-by-line 
transmittance/radiance  code  FASC0D2  (Clough  et  al . ,  1986)  was  exploited  to 
accomplish  this  step  in  formulating  the  physical  retrieval  approach.  FASC0D2 
enables  the  user  to  calculate  high  resolution,  spectrally  dependent 
transmi ttances  and  radiances  by  defining  the  relevant  thermodynamic  and 
compositional  properties  of  the  specific  1 ine - of - s ight  desired.  Thus, 
appropriate  data  on  the  SCRIBE  viewing  geometry  and  first  guesses  for 
temperature  and  emitting  gas  concentrations  could  be  used  to  initialize  the 
retrieval  procedure. 

A  variety  of  minor  changes  were  required  to  tailor  the  application  of  the 
FASC0D2  model  to  our  specific  problem.  FASC0D2 ' s  internal  layering  algorithm 
produced  a  set  of  layers  for  the  vertical  integration  of  radiance  which  was 
inappropriate  for  the  observation  geometry  and  spectral  region  of  interest. 
Selected  layers  were  too  thick  resulting  in  a  suboptimal  vertical 
quadrature.  To  better  accommodate  radiance  calculations  in  the  strongly 
absorbing  spectral  regions  of  interest,  we  apportioned  the  quadrature  levels 
to  be  most  closely  spaced  near  the  float  level  of  30.5  km,  as  follows  for 


downward  viewing:  0(1)28,  28.5,  29,  29.4,  29.7,  30,  30.2,  30.4,  30.5  km;  i.e  , 
36  levels  or  35  layers  in  all,  which  FASC0D2  easily  accommodated. 

We  used  Chahine's  original  form  of  the  relaxation  equation  which,  for  the 
nth  interation,  relates  the  black-body  radiance  B^n^(pj)  for  the  levels  p- , 
corresponding  to  frequency  i/j ,  at  which  temperatures  are  retrieved,  to  the 
previous  blackbody  radiance  B^n'^(pj)  and  to  the  measured  and  calculated 
radiances,  R*  and  R^n’^,  respectively: 

B(n) (p . )  -  B(n"1)(p)R*/R^n"1)  a) 

J  J  J  J 

The  process  is  continued  as  long  as  the  calculated  and  observed  radiances 
continue  to  converge  or  until  their  differences  are  approximately  equal  to  the 
instrumental  and  computational  noise. 

The  15  /jm  band  is  the  appropriate  band  within  the  SCRIBE  spectral  limits 
for  atmospheric  temperature  sensing  (Isaacs  et  al . ,  1986).  Inherent  in  the 
full  resolution  SCRIBE  data  is  considerable  information  pertaining  to 
atmospheric  temperature  structure.  Practical  considerations,  however,  suggest 
the  use  of  a  small  number  of  finite  bandpass  channels  instead  of  the  full 
resolution  spectra.  This  approach  makes  the  retrieval  computationally 
tractable  and  improves  s ignal - to -noise  characteristics  of  the  full  resolution 
data.  There  is  a  trade-off,  however,  between  bandpass  weighting  of  the  full 
resolution  data  and  loss  of  atmospheric  information  inherent  in  the  spectral 
variation  of  radiance.  Therefore,  the  channel  set  selected  has  to  be  optimized 
with  respect  to  these  considerations. 

Such  a  channel  set  for  temperature  sounding  has  been  proposed  in  the 
engineering  feasibility  design  stage  of  the  Advanced  Moisture  and  Temperature 
Sounder  CAMTS)  ,  a  concept  for  future  infrared  satellite  remote  sensing  of  the 
atmosphere.  We  have  adopted  these  preoptimized  channels  in  this  study  of 
temperature  retrieval  from  SCRIBE.  The  AMTS  15  /jm  atmospheric  channels,  are 


given  in  Table  1.  from  Chahine  et  al.  (1984).  These  center  frequencies  were 
used  for  the  SCRIBE  retrieval  and  sensitivity  studies.  SCRIBE  spectra  and 
SCRIBE  spectra  simulated  via  FASC0D2  were  spectrally  integrated  to  the  coarser 
AMTS  channel  resolution  using  triangular  slit  functions  of  0.5  cm'^  full  width 
at  half  height.  Also  given  in  Table  1  are  the  nominal  peaks  of  weighting 
functions  evaluated  for  these  channels  assuming  a  standard  atmosphere  and 
sensor  viewing  nadir  from  space. 


Table  1.  AMTS-based  channel  frequencies  for  SCRIBE  temperature  retrievals. 


Channel 

Number 

Frequency 

(cm-1) 

Halfwidth 

(cm-1) 

Weighting  function 

Peak  (mb) 

1 

606.95 

0.5 

850 

2 

623.20 

0.5 

700 

3 

627.80 

0.5 

400 

4 

634.30 

0.5 

270 

5 

646.60 

0.5 

180 

6 

654.35 

0.5 

90 

7 

665.55 

70 

8 

666.85 

30 

9 

668.15 

0.5 

3 

10 

669.45 

20 

11 

875.00 

900 

The  overall  retrieval  scheme  developed  for  SCRIBE  thus  consists  of  the 

following  steps: 

(a)  Full  resolution  SCRIBE  spectra  are  degraded  with  a  triangular  slit 
function  to  give  channel  radiances  at  0.5  cm-1  resolution  for  the 
sounding  frequencies  in  Table  1;  these  are  the  resampled  data 

(b)  First  guess  path  temperature  and  constituent  profiles  are  assumed  based 
on  a  priori  information  or  climatology  and  FASC0D2  is  run  (with  the 
vertical  layering  modified  for  SCRIBE  runs)  to  calculate  a  full 
resolution  synthetic  spectrum  for  the  wavenumber  domains  of  interest 
(those  included  in  the  channels  of  Table  1. 

(c)  The  full  resolution  spectrum  is  degraded  to  give  0.5  cm-1  channel 
radiances . 


(d)  The  svnthetic  radiances  thus  calculated  are  compared  to  the  degraded 
SCRIBE  data  produced  in  step  (a) . 


(e)  The  temperature  profile  is  adjusted  using  Eqn.  (1). 

(f)  The  process  is  continued  as  long  as  the  calculated  and  observed  radiances 
continue  to  converge  or  until  their  differences  are  approximately  equal 
to  the  instrumental  and  computational  noise. 

The  success  of  the  resulting  retrieval  depends  on  the  levels  of  noise  in 
the  data,  the  quality  of  the  first  guess  and  the  ability  of  the  selected 
channel  set  to  provide  information  on  temperature  variations  along  the 
observational  path.  Since  FASC0D2  is  used  iteratively  and  interactively,  a 
good  first  guess  is  important.  The  efficiency  of  the  retrieval  approach  could 
be  enhanced  by  developing  a  "rapid"  algorithm  to  carry  out  the  path 
transmittance  calculations  required  for  the  forward  problem.  A  rapid 
algorithm  is  essentially  a  band  model  specifically  tuned  for  each  of  the 
desired  synthetic  sensor  channels  (e.g.  Susskind  et  al.,  1982).  It  is  based 
on  an  ensemble  of  calculations  performed  with  the  exact  line-by-line  (LBL) 
code,  i.e.,  FASC0D2 .  This  task  was  not  undertaken  in  Phase  I  due  to  the 
uncertainties  in  treatment  of  line  coupling  effects  at  this  time. 

3  .  2  Retrieval  results 

Two  temperature  retrievals  were  carried  out.  The  first,  to  test  the 
system,  used  the  actual  radiosonde  sounding  to  compute  synthetic  "measured" 
radiances,  and  each  channel  was  associated  with  the  pressure  at  which  the 
channel  brightness  temperature  is  equal  to  the  corresponding  temperature  for  a 
standard  tropical  atmosphere.  An  isothermal  first  guess  was  used,  so  that 
Bj  Rj^n  ^  -  constant  in  Eq .  (1),  and  the  first  iteration  should 

already  give  a  reasonably  good  solution  if  the  actual  sounding  is  close  to 
that  for  the  standard  tropical  atmosphere , which  it  is  in  this  case.  Figure  5 
shows  that  the  solution  has  already  converged  to  the  correct  sounding  after 
only  two  iterations 

The  second  retrieval  exercise  was  carried  out  using  actual  SCRIBE 
measured  radiances.  Data  from  the  nadir  scans  illustrated  in  Figure  1  were 
selected  for  this  test  of  the  retrieval  algorithm.  To  reconcile  apparent 
calibration  problems  in  the  SCRIBE  data,  the  actual  spectra  were  scaled  to 
give  agreement  with  the  blackbody  radiai.ee  corresponding  to  ambient 
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Figure  5  SCRIBE  retrieval  results,  pressure-channel  association  flora  tropical 
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Temperature  at  the  peak  of  the  15  um  Q-branch.  Ambient  temperatures  were 
obtained  from  the  local  radiosonde  data  obtained  in  support  of  the  SCRIBE 
measurements.  Although  the  sounding  was  obtained  a  few  hours  earlier,  there 
should  be  little  diurnal  variation  of  temperature  at  the  levels  observed  at 
the  peak  of  the  Q-branch,  i.e.  essentially  float  altitude. 

In  order  to  adjust  the  guess  profile  using  the  bandpass  weighted  data  and 
simulations  using  Eq.  (1),  it  is  necessary  to  establish  a  one-to-one 
correspondence  between  channel  residuals  and  required  adjustments  at  a 
specified  pressure  level.  The  channels  in  this  case  were  associated  with  the 
pressures  at  which  their  respective  brightness  temperatures  were  equal  to  the 
corresponding  temperatures  of  a  U.S.  Standard  Atmosphere.  This  atmospheric 
model  was  also  employed  as  a  first  guess  of  the  temperature  profile  to 
initiate  the  retrieval  process.  Results  are  shown  in  Figure  6.  Using  the 
radiosonde  as  a  standard  for  comparison  (the  curve  labelled  "SCRIBE"), 
convergence  is  apparent  in  the  general  ranges  between  10  and  30  mb,  200  and 
400  mb,  and  in  the  vicinity  of  the  surface.  Figure  6  also  shows  that  the 
results  were  quite  poor  in  the  lower  stratosphere,  probably  because 
and  R^n  ^  in  Eq .  (1)  are  still  equal,  though  no  longer  constant.  This 
constrains  any  further  relaxation  of  the  lower  stratosphere  temperatures.  One 
further  suspects  that  there  may  be  little  information  on  temperature  variation 
at  these  levels  available  from  the  set  of  selected  channels. 

It  should  be  noted  that  considerable  improvement  can  be  obtained  by 
extrapolating  the  lower  stratospheric  lapse  rate  downwards  and  the  upper 
tropospheric  lapse  rate  upward,  and  taking  the  intersection  as  the  tropopause. 
But  clearly  better  associations  between  channels  and  sounding  levels  are 
required.  We  have  therefore  devoted  the  rest  of  our  remote  sensing  effort  to 
sensitivity  studies  designed  specifically  for  the  SCRIBE  data.  These 
sensitivity  studies  and  their  results  are  described  in  the  following  section. 


3 . 3  Sensitivity  Studies 

The  sensitivity  studies  were  carried  out  in  the  nadir  viewing  mode.  For 
each  set  of  sensitivity  curves,  synthetic  upwelling  SCRIBE  radiances  were 
calculated  for  each  channel  with  the  FASC0D2  forward  program.  The  calcula- 
t i ons  were  done  for  either  a  U.S.  Standard  or  Standard  Tropical  Atmosphere, 
without  perturbation  and  for  each  of  a  set  of  perturbed  atmospheres  obtained 
by  increasing  the  temperature  by  5C  consecutively  over  a  1  km  equivalent  layer 
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Figure  6  SCRIBE  retrieval  results,  pressure -channel  association  from  US 
Standard  Atmosphere : 


thickness  centered  at  the  16  levels  1(2)25(2.5)32.5  km,  then  subtracting  the 
unperturbed  from  the  perturbed  brightness  temperature.  These  brightness 
temperature  differences  are  plotted  as  a  linear  function  of  height  in  Figures 
7  to  16.  They  can  be  used  to  assign  effective  sounding  heights  to  the  various 
channels.  The  sets  of  sensitivity  curves  are  remarkably  similar  to  each 
ocher,  regardless  of  model  atmosphere  used,  which  indicates  that  the  same  set 
of  correspondences  between  channel  and  height  can  be  used  for  all  soundings. 
The  differences  are  discussed  below. 

Figures  7  and  8  show  the  sensitivity  curves  for  the  U.S.  Standard  Atmos¬ 
phere  calculated  from  the  line  parameters  on  the  1982  HITRAN  (Rothman  et  al . , 
1983)  tape,  for  the  troposphere  and  stratosphere,  respectively.  These  can  he 
compared  with  the  sensitivity  curves  calculated  using  the  latest  (1985)  HITRAN 
data  shown  in  Figures  9  and  10.  Line  coupling  was  ignored  in  calculating 
these  sets  of  sensitivity  curves.  The  effect  of  updating  the  line  parameter 
compilation  is  only  evident  for  the  606.95  cm"'''  channel,  which  is  least 
dominated  by  the  lines  of  the  fundamental  transition.  The  606.95  cm'^  channel 
is  redundant  with  the  623.2  cm~  ^  channel;  otherwise  the  tropospheric  channels 
are  nicely  spaced.  The  stratospheric  channels,  on  the  other  hand,  overlap 
markedly.  Further  examining  Figures  7  and  8  it  is  notable  that  there  is  a 
minimum  in  sensitivity  to  temperature  variations  in  the  range  between  about  12 
and  20  km;  the  single  channel  at  654  cm-1  provides  some  information.  This 
corresponds  to  the  pressure  region  between  60  and  200  mb  where  convergence 
problems  were  evident  in  the  retrieval  exercise  described  in  the  previous 
section.  Clearly  further  work  is  required  on  channel  optimization  for  the 
temperature  profile  retrieval  problem. 

Figures  11  and  12  were  calculated  like  Figures  9  and  10,  except  that  it 
was  done  for  the  Standard  Tropical  (M-l  in  FASC0D2)  rather  than  U.S.  Standard 
Atmosphere.  The  tropical  sensitivity  curves  are  indistiguishable  from  the 
standard  curves  in  the  stratosphere,  so  it  shows  the  same  marked  overlapping. 
Limb  scanning  would  be  much  more  appropriate  for  sounding  the  stratosphere. 

The  tropical  sensitivity  curves  for  the  troposphere  are  quite  different  from 
the  standard  curves,  particularly  in  the  region  between  tropical  and  standard 
tropopause,  but  the  heights  of  the  sensitivity  curve  maxima  are  almost 
unchanged.  This  suggests  that  a  universal  association  of  channel  frequence  i 
with  height  z  could  be  chosen  somewhat  as  shown  in  Table  2  on  page  22. 
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Figure  8  SCRIBE  instrument  sensitivity  results:  665-669  cm channels,  old 
line  data,  US  standard  atmosphere,  no  line  coupling. 
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Figure  9  SCRIBE  instrument  sensitivity  results:  606-654  cm'1  channels,  new 
line  data,  US  standard  atmosphere,  no  line  coupling. 
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Figure  11  SCRIBE  instrument  sensitivity  results:  606-654  cm'^  channels 
line  data,  tropical  atmosphere,  no  line  coupling. 
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Figure  12  SCRIBE  instrument  sensitivity  results:  665-669  cm'1  channels 
line  data,  tropical  atmosphere,  no  line  coupling. 


Table  2.  Association  of  channel  frequency  with  altitude  for  retrieval 
application . 

Channel  11  234  5  6  78  9  10 

i/ ( cm’  1 )  875  623.2  627.8  634.3  646.6  654.4  665.6  666.8  668.2  669.4 

z  (km)  13  7  11  14  17  (20)  27.5  29  30 


Sensitivity  curves  were  also  calculated  to  investigate  the  effect  of  line 
coupling  in  the  15  Q-branch.  These  are  presented  in  the  following  section. 

4.  Investigation  of  Line  Coupling  Phenomena  with  Scribe  Data 


4 . 1  Background 


Line  coupling  is  the  distortion  of  the  line  shape  as  a  result  of  the 
rapid  rotational  relaxation  of  closely  spaced  transitions,  such  as  in  Q- 
branches  of  perpendicular  molecules  such  as  CO2 ,  and  perhaps  even  more  so  for 
molecules  with  larger  moments  of  inertia  and  therefore  more  closely  spaced 
lines.  A  discussion  of  this  effect  was  given  together  with  a  bibliography,  by 
M.L.  Hoke  and  S.A.  Clough  of  AFGL  and  W.  Lafferty  and  B.W.  Olson  of  NBS 
(National  Bureau  of  Standards)  at  the  5-6  June  1986  Annual  Review  Conference 
on  Atmospheric  Transmission  Models,  in  a  paper  entitled  "Line  Coupling  in 
Carbon  Dioxide".  Drs .  Hoke  and  Clough  expanded  the  absorption  coefficient  to 
first  order  in  pressure  p,  and  compared  calculated  spectra  of  the  CO2  15  jim  Q- 
branch  with  laboratory  measurements.  Their  line  shape  for  a  line  centered  at 
was  of  the  form 


V"  S  a  +  y  (i/-u  ) 

\  if  if  if  if 

k(v,is  )  =  )  p - - 

if  /  w  2  2  2 

L. -  o  p  +  (u-u  ) 

if  if 

if 


This  expression  reduces  to  the  Lorentz  shape  when  the  coupling  coefficients  Y 
are  set  to  zero.  The  Y's  used  in  the  laboratory  comparison  were  calculated 
for  the  laboratory  temperature  of  296  K,  and  are  included  in  the  latest  HITRAN 


tape  (Rothman,  1986).  The  introduction  of  line  coupling  almost  eliminated  a 
rather  sizable  discrepancy  between  observed  spectra  and  spectra  calculated 
with  the  use  of  a  Lorentz  line  shape. 


CO2  line  coupling  near  15  has  nor  been  investigated  at  other  than 

laboratory  room  temperature.  The  SCRIBE  data,  on  the  other  hand,  provide  the 
opportunity  to  check  calculations  of  the  temperature  dependence  of  line 
coupling,  and  its  effect  on  atmospheric  t ransmi r r ance  and  radiances.  This  is 
especially  important  in  the  main  Q-branch  region,  which  is  useful  for 
retrieval  of  stratospheric  temperatures. 

Dr.  Hoke  kindly  agreed  to  calculate  Y  values  for  the  temperatures  210, 
230,  260,  296,  330  ,  so  that  we  could  investigate  the  effect  of  temperature 
dependence  variations  on  radiances  and  on  the  sensitivity  curves. 

4 • 2  Effect  of  Line  Coupling  on  Sensitivity  Curves 

Sensitivity  curves  were  calculated  with  line  coupling  included  by  using 
the  line  coupling  option  of  FASC0D2.  As  currently  formulated  line  coupling 
coefficients  are  stored  in  the  space  allocated  in  the  line  parameter 
compilation  for  pressure  shifts.  When  non-zero  values  are  encountered  for 
this  parameter  the  FASC0D2  algorithm  calculates  the  appropriate  correction  to 
line  shape  using  Eq.  (2)  above.  To  accomodate  Hoke's  temperature  dependent 
line  coupling  coefficients,  the  line  parameter  compilation  was  edited  to 
substitute  the  more  recent  values  at,  for  example,  230K.  Since  only  the  Y's 
for  the  fundamental  Q-branch  were  used,  the  major  effects  could  be  expected  to 
be  found  for  the  channels  sounding  just  below  SCRIBE  float  altitude  where  the 
temperature  230K  was  representative. 

It  is  important  to  note  that  this  application  of  FASC0D2  uses  the  line 
coupling  parameters  for  the  chosen  temperature  for  all  layers,  i.e.  at  all 
selected  mean  layer  temperatures  within  the  simulation.  The  calculation  is 
therefore  isothermal  and  does  nor  fully  exploit  the  capability  to  calculate 
temperature  dependent  line  coupling  coefficients.  This  can  be  accomp 1 i s hed , 
bit  only  through  an  appropriate  modification  of  the  FASC0D2  algorithm  to 
either  accept  temperature  dependent  line  parameter  data  (i.e.  line  parameters 
f  ;r  each  layer)  or  by  fitting  a  simple  functional  form  to  or  sealing  the  line 
coupling  coefficient  temperature  dependence  so  that  it  can  be  internally 
calculated  as  is  currently  done  for  line  strength,  for  example. 


The  tropospheric  and  stratospheric  results  with  line  coupling, 
respectively,  are  shown  in  Figures  13  and  14  for  the  US.  standard  atmosphere 
and  in  Figures  13  and  16,  respectively,  for  the  standard  tropical  (M-l) 
atmosphere.  When  these  sensitivity  curves  are  compared  with  Figures  9  through 
12  it  is  seen  that  potential  effects  of  line  coupling  are  evident  on  the 
sensitivity  curves  for  the  666.8  cm"''-  channel,  but  that  the  sensitivity  curves 
for  the  other  channels  are  relatively  unaffected.  To  see  why  this  is  so  for 
the  stratospheric  channels  we  must  look  at  the  spectrally  detailed  synthetic 
radiances  shown  in  the  figures  to  follow.  The  apparent  insensitivity  of  the 
tropospheric  channel  radiances  to  line  coupling  is  to  be  expected  because  we 
ignored  all  Q-branches  except  for  the  fundamental.  The  entire  15  ^m  band 
region  is  strewn  with  Q-branches  separated  from  each  other  by  accidental  Fermi 
resonances.  Thus  the  effect  found  for  the  666.8  cm'^  channel  may  also  appear 
in  tropospheric  or  tropopause  channels  such  as  that  at  646.6  cm  ^ ,  which  is 
very  close  to  two  Q-branches.  This  emphasizes  the  importance  of  line  coupling 
to  remote  sensing  and  the  necessity  for  extending  the  investigation  of  these 
sensitivities  to  other  related  Q-branches. 

It  should  also  be  noted  that  even  for  the  666.8  cm~  ^  channel,  the  height 
of  the  peak  of  the  sensitivity  curve  is  unaffected  by  line  coupling. 


Comparison  of  SCRIBE  Data  to  Synthetic  Calculations 


Detailed  synthetic  SCRIBE  spectra  were  calculated  using  the  modified 
FASC0D2  algorithm  from  660  to  670  cm"V  Line  coupling  effects  were  not 
included  in  the  calculation  due  to  the  non  isothermal  nature  of  the  emission 
path.  Nadir  brightness  temperatures  and  radiances  are  shown  in  Figures  17  and 
18,  respectively.  These  can  be  compared  to  the  corresponding  observed  SCRIBE 
nadir  viewing  spectra  in  Figures  19  and  20  (e.g  Kelvins  and  radiance, 
respectively.)  Here  there  is  an  obvious  calibration  problem,  since  the 
brightness  temperature  measured  at  the  Q-branch  maximum  is  close  to  250K 
rather  than  the  ambient  temperature  of  230K.  When  the  measured  radiances  are 
shifted  downward  on  a  logarithmic  radiance  scale  to  peak  at  230K  at  the 
maxima,  and  a  window  correction  of  five  percent  of  ambient  temperature  is 
added,  the  resulting  spectrum  is  reasonably  similar  to  that  for  the  synthetic 
spectra.  This  assumes  that  the  synthetic  spectra  are  first  degraded  in 
resolution  to  that  of  the  SCRIBE  instrument. 
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Figure  13  SCRIBE  instrument  sensitivity  results.  606-669  cm’1  channels,  ole 
line  data,  US  standard  atmosphere,  with  line  coupling  calculated  f 
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Figure  16  SCRIBE  instrument  sensitivity  results  665-669  cm'^  channels,  new 
line  data,  tropical  atmosphere,  with  line  coupling  calculated  for 
?  30K . 
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To  investigate  potential  effects  of  line  coupling  on  the  spectra  in  the 
immediate  vicinity  of  the  Q  branch,  synthetic  spectra  were  evaluated  with  and 
without  line  coupling  effects  in  the  reduced  wavenumber  domain  from  666.5  to 
668.0  cm'V  This  spectrum  (without  line  coupling)  is  illustrated  in  Figure 
21.  Figure  22  shows  the  same  result,  however  including  line  coupling  effects 
based  on  Hoke's  coupling  coefficients  for  210K.  The  reduced  wavenumber  domain 
was  chosen  to  isolate  a  spectral  region  whose  emission  could  be  characterized 
by  assuming  a  quasi- isothermal  layer  at  210K  to  facilitate  the  line  coupling 
calculation  The  simulations  with  and  without  line  coupling  (Figures  22  and 
21,  respectively)  can  be  compared  to  the  observed  SCRIBE  spectrum  (Figure 
23).  The  reduced  radiance  from  666.5  to  667.3  cm"'*'  can  be  attributed  to  a 
decreased  absorption  coefficient  induced  by  line  coupling  in  the  part  of  the 
stratosphere  below  the  SCRIBE  level  where  the  temperature  increases  with 
height,  i.e.,  SCRIBE  "sees"  less  warm  and  more  cold  air  because  of  line 
coupling.  These  results  help  explain  the  effects  of  line  coupling  on  the 
sensitivity  curves.  The  666.8  cm’*  channel  is  centrally  located  in  the 
spectral  region  where  there  is  maximum  weakening  of  the  absorption  coefficient 
from  line  coupling  and  therefore  should  show,  as  it  does,  greatest  decrease  in 
sensitivity.  The  channels  at  668.2  and  669.4  cm’*  are  so  opaque  that  even  a 
large  relative  change  in  absorption  coefficient  should  not  affect  the  opacity. 


Figures  24  to  25  present  detailed  synthetic  spectra  as  in  Figures  17  and 
18.  but  for  limb  viewing  at  a  zenith  angle  of  93.7°.  The  corresponding 
measured  SCRIBE  spectra,  shown  in  Figures  26  and  27,  do  not  seem  to  have  the 
residual  calibration  problems  that  the  nadir  spectra  do.  Again,  the 
calculated  line  coupling  effects  are  evidenced  by  a  reduced  radiance.  This 
can  be  seen  by  comparing  Figure  28  and  29  (no  line  coupling  and  line  coupling, 
respectively).  The  agreement  between  the  observed  spectrum  (Figure  30)  and 
the  calculation  is  improved  by  the  incorporation  of  the  line  coupling 
effect.  In  this  case,  that  is  to  be  expected  because  the  viewing  tangent 
height  is  about  17  km,  so  that  most  of  the  trajectory  is  through  an 
atmospheric  path  at  heights  where  the  temperature  is  less  than  210K.  In  the 
opaque  spectral  regions,  where  the  brightness  temperatures  are  higher,  the 
tangent  height  temperature  is  irrelevant,  as  are  line  coupling  effects. 

For  limb  scanning,  th-_  agreement  between  calculated  and  measured  spectra 
is  at  least  as  good  as  for  nadir  viewing.  Again  the  effects  are  potentially 
very  important,  and  again  the  temperature  dependence  of  these  effects  are  even 
more  so . 


TEMPERATURE 


WAVENUMBER  (CM- 1 ) 


Figure  21  FASCODE  simulation  of  SCRIBE  nadir  data  for  the  666.5-668  cm'1  range 
(no  line  coupling).  (in  deg  K) 


RADIANCE  WATTS  /  (CM*«2*STER*CM- 1 ) 
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15-12 


WAVELENGTH  (MICROMETER) 

15.08  15.04  15.00  14.96 


WAVENUMBER  (CM- 1 ) 


Figure  25  FASCODE  simulation  of  SCRIBE  limb  data  for  the  660-670  cm 
no  line  coupling).  (radiance) 
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RADIANCE  WATTS  /  (CM**2«STER*CM- 1 ) 


SCRIBE  DATA  (93.7  DEGREES) 


(in  dog  K ) 


b  Conclusions  and  Recommendations 

!  S  ■  imma  r  v 

WV  have  developed  a  FASCODE  routine  for  calculating  synthetic  radiances 
:  i.at  would  be  received  at  SCRIBE  level  and  viewing  angle,  and  used  this 
so:  i:it'  to  conduct  temperature  retrieval  studies  and  to  investigate  the  effect 
•t  line  coupling  on  radiances. 

Physical  retrievals  obtained  with  the  use  of  the  AMTS  15  /rm  channels  were 
tound  to  he  sensitive  to  the  first  guess  to  the  solution.  A  series  of 
setts i t ivit v  studies  were  carried  out  and  resulted  in  finding  a  nicely  spaced 
universal  set  of  heights  to  which  the  tropospheric  channels  could  correspond 
i :  t  : e laxat ion- tvpe  retrievals.  We  1 1  -  resolved  stratospheric  temperature 
retrievals  probably  require  limb  scanning. 

bine  coupling  was  found  to  have  strong  effects  on  atmospheric  radiances 
ind  to  be  temperature  sensitive.  It  seems  to  have  the  potential  for  being 
;  i r*:cu 1 ar 1 v  important  in  the  troposphere,  and  may  be  useful  for  remote 
1  u  i-  ion.  of  the  tropopause 

rhe  SCRIBE  spectra  are  ideal  for  these  types  of  study. 

Re  *  ndat  ions 

It  is  stronglv  recommended  that  study  and  interpretation  of  the  SCRIBE 
it*.,  he  continued  and  expanded,  with  emphasis  on  both  retrieval  studies  and 

1  i  coupling  effects. 

More  work  should  he  done  on  calibrating  the  spectra,  and  the  new  1986 
;.i  : i  ,i  should  he  calibrated,  studied  and  used.  It  is  clear  from  our 
•  >trp.»r  i  son  between  synthetic  and  actual  spectra  that  there  are  discrepancies 
which  require  further  investigation.  More  SCRIBE  spectra  should  be  taken,  and 
ext mded  to  other  spectral  regions.  There  is  much  to  be  obtained  from  them. 

: he  retrieval  studies  undertaken  in  this  initial  phase  of  our  research 
provide  the  baseline  capabilities  for  exploiting  the  extensive  information 
content  of  the  SCRIBE  spectra.  The  sensitivity  calculations  have  provided  the 
background  for  selection  of  an  optimal  channel  set  for  temperature  retrieval 
studies  Limited  resources  in  the  initial  effort  preclude  taking  advantage  of 
the  temperature  retrieval  capability  at  this  time  Similar  considerations 


apply  Co  the  development  of  analogous  constituent  retrieval  capabilities  for 
which  the  temperature  retrieval  is  a  necessary  prerequisite.  Approaches  have 
been  developed  for  constituent  channel  selection  and  related  sensitivity 
studies . 

To  significantly  improve  retrieval  efficiency,  rapid  algorithm  retrieval 
routines  should  be  developed  for  the  specified  SCRIBE  sounding  channels.  To 
enhance  the  applicability  of  the  FASC0D2  forward  problem  calculation,  variabl 
temperature  dependence  of  the  line  coupling  coefficients  should  be 
incorporated  in  the  FASCODE  program  and  also  in  the  rapid  algorithm  routines, 
if  possible,  and  programs  should  be  developed  for  obtaining  temperature  and 
constituent  retrievals  from  limb  sounding. 

Examination  of  line  coupling  effects  should  be  extended  to  other  CO2  Q- 
branches  in  the  15  ^m  complex  and  to  Q-brances  of  freons ,  C  NOj,  0-j ,  etc.  Th 
effect  of  line  coupling  on  remote  sensing  of  temperatue  and  constituent 
concentrations  should  be  investigated  further  taking  advantage  of  the  FASC0D2 
based  retrieval  capability  developed  here  and  its  future  enhancements. 
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